A concerted neuron-astrocyte program declines in ageing and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Concerted neuron-astrocyte gene expression declines in aging and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a striking relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA-seq to analyze the prefrontal cortex of 191 human donors ages 22-97 years, including healthy individuals and persons with schizophrenia. Latent-factor analysis of these data revealed that in persons whose cortical neurons more strongly expressed genes for synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the Synaptic Neuron-and-Astrocyte Program (SNAP). In schizophrenia and aging - two conditions that involve declines in cognitive flexibility and plasticity 1,2 - cells had divested from SNAP: astrocytes, glutamatergic (excitatory) neurons, and GABAergic (inhibitory) neurons all reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy persons of similar age, may underlie many aspects of normal human interindividual differences and be an important point of convergence for multiple kinds of pathophysiology.

Natural variation in gene expression and viral susceptibility revealed by neural progenitor cell villages.

Human genome variation contributes to diversity in neurodevelopmental outcomes and vulnerabilities; recognizing the underlying molecular and cellular mechanisms will require scalable approaches. Here, we describe a "cell village" experimental platform we used to analyze genetic, molecular, and phenotypic heterogeneity across neural progenitor cells from 44 human donors cultured in a shared in vitro environment using algorithms (Dropulation and Census-seq) to assign cells and phenotypes to individual donors. Through rapid induction of human stem cell-derived neural progenitor cells, measurements of natural genetic variation, and CRISPR-Cas9 genetic perturbations, we identified a common variant that regulates antiviral IFITM3 expression and explains most inter-individual variation in susceptibility to the Zika virus. We also detected expression QTLs corresponding to GWAS loci for brain traits and discovered novel disease-relevant regulators of progenitor proliferation and differentiation such as CACHD1. This approach provides scalable ways to elucidate the effects of genes and genetic variation on cellular phenotypes.

Efficient generation of lower induced motor neurons by coupling Ngn2 expression with developmental cues.

Human pluripotent stem cells (hPSCs) are a powerful tool for disease modeling of hard-to-access tissues (such as the brain). Current protocols either direct neuronal differentiation with small molecules or use transcription-factor-mediated programming. In this study, we couple overexpression of transcription factor Neurogenin2 (Ngn2) with small molecule patterning to differentiate hPSCs into lower induced motor neurons (liMoNes/liMNs). This approach induces canonical MN markers including MN-specific Hb9/MNX1 in more than 95% of cells. liMNs resemble bona fide hPSC-derived MN, exhibit spontaneous electrical activity, express synaptic markers, and can contact muscle cells in vitro. Pooled, multiplexed single-cell RNA sequencing on 50 hPSC lines reveals reproducible populations of distinct subtypes of cervical and brachial MNs that resemble their in vivo, embryonic counterparts. Combining small molecule patterning with Ngn2 overexpression facilitates high-yield, reproducible production of disease-relevant MN subtypes, which is fundamental in propelling our knowledge of MN biology and its disruption in disease.

Molecular convergence between Down syndrome and fragile X syndrome identified using human pluripotent stem cell models.

Down syndrome (DS), driven by an extra copy of chromosome 21 (HSA21), and fragile X syndrome (FXS), driven by loss of the RNA-binding protein FMRP, are two common genetic causes of intellectual disability and autism. Based upon the number of DS-implicated transcripts bound by FMRP, we hypothesize that DS and FXS may share underlying mechanisms. Comparing DS and FXS human pluripotent stem cell (hPSC) and glutamatergic neuron models, we identify increased protein expression of select targets and overlapping transcriptional perturbations. Moreover, acute upregulation of endogenous FMRP in DS patient cells using CRISPRa is sufficient to significantly reduce expression levels of candidate proteins and reverse 40% of global transcriptional perturbations. These results pinpoint specific molecular perturbations shared between DS and FXS that can be leveraged as a strategy for target prioritization; they also provide evidence for the functional relevance of previous associations between FMRP targets and disease-implicated genes.

Whole-genome analysis of human embryonic stem cells enables rational line selection based on genetic variation.

Despite their widespread use in research, there has not yet been a systematic genomic analysis of human embryonic stem cell (hESC) lines at a single-nucleotide resolution. We therefore performed whole-genome sequencing (WGS) of 143 hESC lines and annotated their single-nucleotide and structural genetic variants. We found that while a substantial fraction of hESC lines contained large deleterious structural variants, finer-scale structural and single-nucleotide variants (SNVs) that are ascertainable only through WGS analyses were present in hESC genomes and human blood-derived genomes at similar frequencies. Moreover, WGS allowed us to identify SNVs associated with cancer and other diseases that could alter cellular phenotypes and compromise the safety of hESC-derived cellular products transplanted into humans. As a resource to enable reproducible hESC research and safer translation, we provide a user-friendly WGS data portal and a data-driven scheme for cell line maintenance and selection.

Publicly Available hiPSC Lines with Extreme Polygenic Risk Scores for Modeling Schizophrenia.

Schizophrenia (SZ) is a common and debilitating psychiatric disorder with limited effective treatment options. Although highly heritable, risk for this polygenic disorder depends on the complex interplay of hundreds of common and rare variants. Translating the growing list of genetic loci significantly associated with disease into medically actionable information remains an important challenge. Thus, establishing platforms with which to validate the impact of risk variants in cell-type-specific and donor-dependent contexts is critical. Towards this, we selected and characterized a collection of 12 human induced pluripotent stem cell (hiPSC) lines derived from control donors with extremely low and high SZ polygenic risk scores (PRS). These hiPSC lines are publicly available at the California Institute for Regenerative Medicine (CIRM). The suitability of these extreme PRS hiPSCs for CRISPR-based isogenic comparisons of neurons and glia was evaluated across 3 independent laboratories, identifying 9 out of 12 meeting our criteria. We report a standardized resource of publicly available hiPSCs on which we hope to perform genome engineering and generate diverse kinds of functional data, with comparisons across studies facilitated by the use of a common set of genetic backgrounds.

The genetic architecture of DNA replication timing in human pluripotent stem cells.

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome's replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) - sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Publisher Correction: Exome sequencing in amyotrophic lateral sclerosis implicates a novel gene, DNAJC7, encoding a heat-shock protein.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing in amyotrophic lateral sclerosis implicates a novel gene, DNAJC7, encoding a heat-shock protein.

To discover novel genes underlying amyotrophic lateral sclerosis (ALS), we aggregated exomes from 3,864 cases and 7,839 ancestry-matched controls. We observed a significant excess of rare protein-truncating variants among ALS cases, and these variants were concentrated in constrained genes. Through gene level analyses, we replicated known ALS genes including SOD1, NEK1 and FUS. We also observed multiple distinct protein-truncating variants in a highly constrained gene, DNAJC7. The signal in DNAJC7 exceeded genome-wide significance, and immunoblotting assays showed depletion of DNAJC7 protein in fibroblasts in a patient with ALS carrying the p.Arg156Ter variant. DNAJC7 encodes a member of the heat-shock protein family, HSP40, which, along with HSP70 proteins, facilitates protein homeostasis, including folding of newly synthesized polypeptides and clearance of degraded proteins. When these processes are not regulated, misfolding and accumulation of aberrant proteins can occur and lead to protein aggregation, which is a pathological hallmark of neurodegeneration. Our results highlight DNAJC7 as a novel gene for ALS.

Dysregulated protocadherin-pathway activity as an intrinsic defect in induced pluripotent stem cell-derived cortical interneurons from subjects with schizophrenia.

We generated cortical interneurons (cINs) from induced pluripotent stem cells derived from 14 healthy controls and 14 subjects with schizophrenia. Both healthy control cINs and schizophrenia cINs were authentic, fired spontaneously, received functional excitatory inputs from host neurons, and induced GABA-mediated inhibition in host neurons in vivo. However, schizophrenia cINs had dysregulated expression of protocadherin genes, which lie within documented schizophrenia loci. Mice lacking protocadherin-α showed defective arborization and synaptic density of prefrontal cortex cINs and behavioral abnormalities. Schizophrenia cINs similarly showed defects in synaptic density and arborization that were reversed by inhibitors of protein kinase C, a downstream kinase in the protocadherin pathway. These findings reveal an intrinsic abnormality in schizophrenia cINs in the absence of any circuit-driven pathology. They also demonstrate the utility of homogenous and functional populations of a relevant neuronal subtype for probing pathogenesis mechanisms during development.

Convergence of independent DISC1 mutations on impaired neurite growth via decreased UNC5D expression.

The identification of convergent phenotypes in different models of psychiatric illness highlights robust phenotypes that are more likely to be implicated in disease pathophysiology. Here, we utilize human iPSCs harboring distinct mutations in DISC1 that have been found in families with major mental illness. One mutation was engineered to mimic the consequences on DISC1 protein of a balanced translocation linked to mental illness in a Scottish pedigree; the other mutation was identified in an American pedigree with a high incidence of mental illness. Directed differentiation of these iPSCs using NGN2 expression shows rapid conversion to a homogenous population of mature excitatory neurons. Both DISC1 mutations result in reduced DISC1 protein expression, and show subtle effects on certain presynaptic proteins. In addition, RNA sequencing and qPCR showed decreased expression of UNC5D, DPP10, PCDHA6, and ZNF506 in neurons with both DISC1 mutations. Longitudinal analysis of neurite outgrowth revealed decreased neurite outgrowth in neurons with each DISC1 mutation, which was mimicked by UNC5D knockdown and rescued by transient upregulation of endogenous UNC5D. This study shows a narrow range of convergent phenotypes of two mutations found in families with major mental illness, and implicates dysregulated netrin signaling in DISC1 biology.

Oligodendrocyte differentiation of induced pluripotent stem cells derived from subjects with schizophrenias implicate abnormalities in development.

Abnormalities of brain connectivity and signal transduction are consistently observed in individuals with schizophrenias (SZ). Underlying these anomalies, convergent in vivo, post mortem, and genomic evidence suggest abnormal oligodendrocyte (OL) development and function and lower myelination in SZ. Our primary hypothesis was that there would be abnormalities in the number of induced pluripotent stem (iPS) cell-derived OLs from subjects with SZ. Our secondary hypothesis was that these in vitro abnormalities would correlate with measures of white matter (WM) integrity and myelination in the same subjects in vivo, estimated from magnetic resonance imaging. Six healthy control (HC) and six SZ iPS cell lines, derived from skin fibroblasts from well-characterized subjects, were differentiated into OLs. FACS analysis of the oligodendrocyte-specific surface, glycoprotein O4, was performed at three time points of development (days 65, 75, and 85) to quantify the number of late oligodendrocyte progenitor cells (OPCs) and OLs in each line. Significantly fewer O4-positive cells developed from SZ versus HC lines (95% CI 1.0: 8.6, F1,10 = 8.06, p = 0.02). The difference was greater when corrected for age (95% CI 5.4:10.4, F1,8 = 53.6, p < 0.001). A correlation between myelin content in WM in vivo, estimated by magnetization transfer ratio (MTR) and number of O4-positive cells in vitro was also observed across all time points (F1,9 = 4.3, p = 0.07), reaching significance for mature OLs at day 85 in culture (r = 0.70, p < 0.02). Low production of OPCs may be a contributing mechanism underlying WM reduction in SZ.

The C9orf72-interacting protein Smcr8 is a negative regulator of autoimmunity and lysosomal exocytosis.

While a mutation in C9ORF72 is the most common genetic contributor to amyotrophic lateral sclerosis (ALS), much remains to be learned concerning the function of the protein normally encoded at this locus. To elaborate further on functions for C9ORF72, we used quantitative mass spectrometry-based proteomics to identify interacting proteins in motor neurons and found that its long isoform complexes with and stabilizes SMCR8, which further enables interaction with WDR41. To study the organismal and cellular functions for this tripartite complex, we generated Smcr8 loss-of-function mutant mice and found that they developed phenotypes also observed in C9orf72 loss-of-function animals, including autoimmunity. Along with a loss of tolerance for many nervous system autoantigens, we found increased lysosomal exocytosis in Smcr8 mutant macrophages. In addition to elevated surface Lamp1 (lysosome-associated membrane protein 1) expression, we also observed enhanced secretion of lysosomal components-phenotypes that we subsequently observed in C9orf72 loss-of-function macrophages. Overall, our findings demonstrate that C9ORF72 and SMCR8 have interdependent functions in suppressing autoimmunity as well as negatively regulating lysosomal exocytosis-processes of potential importance to ALS.

Combining NGN2 Programming with Developmental Patterning Generates Human Excitatory Neurons with NMDAR-Mediated Synaptic Transmission.

Transcription factor programming of pluripotent stem cells (PSCs) has emerged as an approach to generate human neurons for disease modeling. However, programming schemes produce a variety of cell types, and those neurons that are made often retain an immature phenotype, which limits their utility in modeling neuronal processes, including synaptic transmission. We report that combining NGN2 programming with SMAD and WNT inhibition generates human patterned induced neurons (hpiNs). Single-cell analyses showed that hpiN cultures contained cells along a developmental continuum, ranging from poorly differentiated neuronal progenitors to well-differentiated, excitatory glutamatergic neurons. The most differentiated neurons could be identified using a CAMK2A::GFP reporter gene and exhibited greater functionality, including NMDAR-mediated synaptic transmission. We conclude that utilizing single-cell and reporter gene approaches for selecting successfully programmed cells for study will greatly enhance the utility of hpiNs and other programmed neuronal populations in the modeling of nervous system disorders.

Human pluripotent stem cells recurrently acquire and expand dominant negative P53 mutations.

Human pluripotent stem cells (hPS cells) can self-renew indefinitely, making them an attractive source for regenerative therapies. This expansion potential has been linked with the acquisition of large copy number variants that provide mutated cells with a growth advantage in culture. The nature, extent and functional effects of other acquired genome sequence mutations in cultured hPS cells are not known. Here we sequence the protein-coding genes (exomes) of 140 independent human embryonic stem cell (hES cell) lines, including 26 lines prepared for potential clinical use. We then apply computational strategies for identifying mutations present in a subset of cells in each hES cell line. Although such mosaic mutations were generally rare, we identified five unrelated hES cell lines that carried six mutations in the TP53 gene that encodes the tumour suppressor P53. The TP53 mutations we observed are dominant negative and are the mutations most commonly seen in human cancers. We found that the TP53 mutant allelic fraction increased with passage number under standard culture conditions, suggesting that the P53 mutations confer selective advantage. We then mined published RNA sequencing data from 117 hPS cell lines, and observed another nine TP53 mutations, all resulting in coding changes in the DNA-binding domain of P53. In three lines, the allelic fraction exceeded 50%, suggesting additional selective advantage resulting from the loss of heterozygosity at the TP53 locus. As the acquisition and expansion of cancer-associated mutations in hPS cells may go unnoticed during most applications, we suggest that careful genetic characterization of hPS cells and their differentiated derivatives be carried out before clinical use.

Loss-of-function mutations in the C9ORF72 mouse ortholog cause fatal autoimmune disease.

C9ORF72 mutations are found in a significant fraction of patients suffering from amyotrophic lateral sclerosis and frontotemporal dementia, yet the function of the C9ORF72 gene product remains poorly understood. We show that mice harboring loss-of-function mutations in the ortholog of C9ORF72 develop splenomegaly, neutrophilia, thrombocytopenia, increased expression of inflammatory cytokines, and severe autoimmunity, ultimately leading to a high mortality rate. Transplantation of mutant mouse bone marrow into wild-type recipients was sufficient to recapitulate the phenotypes observed in the mutant animals, including autoimmunity and premature mortality. Reciprocally, transplantation of wild-type mouse bone marrow into mutant mice improved their phenotype. We conclude that C9ORF72 serves an important function within the hematopoietic system to restrict inflammation and the development of autoimmunity.

Genetic variation in human DNA replication timing.

Genomic DNA replicates in a choreographed temporal order that impacts the distribution of mutations along the genome. We show here that DNA replication timing is shaped by genetic polymorphisms that act in cis upon megabase-scale DNA segments. In genome sequences from proliferating cells, read depth along chromosomes reflected DNA replication activity in those cells. We used this relationship to analyze variation in replication timing among 161 individuals sequenced by the 1000 Genomes Project. Genome-wide association of replication timing with genetic variation identified 16 loci at which inherited alleles associate with replication timing. We call these "replication timing quantitative trait loci" (rtQTLs). rtQTLs involved the differential use of replication origins, exhibited allele-specific effects on replication timing, and associated with gene expression variation at megabase scales. Our results show replication timing to be shaped by genetic polymorphism and identify a means by which inherited polymorphism regulates the mutability of nearby sequences.

Sensory maps in the olfactory cortex defined by long-range viral tracing of single neurons.

Sensory information may be represented in the brain by stereotyped mapping of axonal inputs or by patterning that varies between individuals. In olfaction, a stereotyped map is evident in the first sensory processing centre, the olfactory bulb (OB), where different odours elicit activity in unique combinatorial patterns of spatially invariant glomeruli. Activation of each glomerulus is relayed to higher cortical processing centres by a set of ∼20-50 'homotypic' mitral and tufted (MT) neurons. In the cortex, target neurons integrate information from multiple glomeruli to detect distinct features of chemically diverse odours. How this is accomplished remains unclear, perhaps because the cortical mapping of glomerular information by individual MT neurons has not been described. Here we use new viral tracing and three-dimensional brain reconstruction methods to compare the cortical projections of defined sets of MT neurons. We show that the gross-scale organization of the OB is preserved in the patterns of axonal projections to one processing centre yet reordered in another, suggesting that distinct coding strategies may operate in different targets. However, at the level of individual neurons neither glomerular order nor stereotypy is preserved in either region. Rather, homotypic MT neurons from the same glomerulus innervate broad regions that differ between individuals. Strikingly, even in the same animal, MT neurons exhibit extensive diversity in wiring; axons of homotypic MT pairs diverge from each other, emit primary branches at distinct locations and 70-90% of branches of homotypic and heterotypic pairs are non-overlapping. This pronounced reorganization of sensory maps in the cortex offers an anatomic substrate for expanded combinatorial integration of information from spatially distinct glomeruli and predicts an unanticipated role for diversification of otherwise similar output neurons.